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A new microscopic model of the charge Kondo effect is proposed, and is studied on the basis 
of a numerical renormalization group calculations, ft is shown that the charge Kondo effect is 
caused by a pair hopping interaction between the conduction band and the localized band. It is 
found that the formation of a charge Kondo- Yosida singlet due to its mechanism corresponds to 
the appearance of valence skipping phenomenon, and that the charge Kondo- Yosida singlet and 
the spin Kondo- Yosida singlet state coexist for a tuned parameter set, because the charge and 
spin degrees of freedom are separated in this model. Then, it is also found that the Sommerfeld 
coefficient is enhanced by the competition between these two singlet states. 

KEYWORDS: Tl doped PbTe, Charge Kondo Effect, Pair Hopping, Negative-U Mechanism, Numerical 
renormalization group 



It is well known that the formal valence of a Tl ion is 
1+ or 3 +, 2+ being skipped, and that of Sn and Pb is 
2+ and 4+, 3+ being skipped. The same phenomenon 
also appears for As, Bi, and so on. Namely, the formal 
number of ns electrons of these ions is either (ns ) or 
(ns 2 ), while an (ns 1 ) state does not appear, where n is 
the principal quantum number of an electronic state in 
centro-symmetric ions. This phenomenon is called "va- 
lence skipping" or "inert pair effect" . 

To explain this phenomenon, many mechanisms have 
been proposed so far: An intra-atomic mechanism, an 
inter-atomic mechanism, a strong electron-phonon cou- 
pling mechanism, and so on. 1-7 It has been found that 
the phenomenon can be described by the attractive on- 
site interaction (i.e., negative-U model). Thus, on the 
basis of the negative-U model, the electronic states of 
compounds including the valence skipping elements have 
been studied, without asking the origin of the negative- 
U. 

Recently, a new type of Kondo effect has been observed 
in Tl doped PbTe. 8 Since Tl impurity has no magnetic 
moment, it has been suggested that the charge degener- 
acy of Tl ion is important for such a Kondo like effect 
to appear. Thus, this effect is called the charge Kondo 
effect. To explain the origin of the charge Kondo effect 
theoretically, the negative-U Anderson model 9 has been 
proposed, and the electronic state has been discussed in 
detail. 10 ' 11 

However, there are few discussions about the charge 
Kondo effect and the valence skipping phenomenon based 
on a microscopic model without assuming the negative- 
U. In addition, although novel physical properties are 
expected by the competition of the charge Kondo effect 
and the spin Kondo effect, it has not been discussed in 
detail. 

In this letter, we propose a microscopic mechanism 
based on electron correlations: It is shown that the charge 
Kondo effect is caused by a pair hopping interaction be- 
tween a conduction band and a localized state, and is 



found that the charge Kondo- Yosida (KY) singlet and 
the spin KY singlet can coexist for a certain set of pa- 
rameters in the microscopic model. 

We propose an effective model Hamiltonian as follows: 



tl — Ho + rlU + Hhybi 

where Ho = H c + Hd + Hdc + H p h with 
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Here, d a and c k(J are annihilation operators of an elec- 
tron on the localized state and the conduction band with 
wave vector k and spin a, and n^a = d\d a . We denote 
the localized state as "d" state. £k and u are the dis- 
persion of conduction band and the chemical potential, 
respectively. A^, Udc, and J v h are the one-body level of 
the localized state, the inter-orbital Coulomb interaction, 
and the pair hopping interaction between the conduc- 
tion band and the localized state, respectively. Ud and 
Vdc are the intra-orbital Coulomb interaction on the lo- 
calized state and the hybridization between the localized 
state and the conduction band. 

To clarify a role of the charge degree of freedom, we 
introduce pseudo-spin representations as follows: 



1 . 

~{n dt + n di - 1), 



(8) 
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Hereafter, quantities denned by eqs. (8) ~ (13) are called 
axial charges. Then, the part T~Lq of model Hamiltonian 
(1) is transformed as 

Ho^H a = H c + U P ot + iidc + H P h + Hd, (14) 
where 



J^(ek - M + t / dc)c ] | C(T Ck ', (15) 

k<r 



H 



pot 



H dc 



kk'(k^k') 

W dc P d I Z c , 



(9) The initial Hamiltonian H$ is given by 
(10) #o = A- 1/2 [Hdc + n p h + Hd\- (22) 

Here we have discarded T-L po t, because H po t oa U dc (l — 
A -1 )— s-0 for A— > 1 according to the discussion in ref. 13 
One can estimate the eigenvalues and eigenstates of 
Hm by the repeated use of the recursion form with keep- 
ing 100 states in each iteration step. We also performed 
the NRG calculations with keeping 300 states, and ver- 
ified that the result was the same as that keeping 100 
states. 

In this letter, we show the entropy of the localized 
state and the fraction of the zero-, single-, and double 
occupancy on the localized state (Dq d ), where (DQ d ) 
takes Qd = (for d° state), Q d = 1 (for d 1 state), and 
Q d = 2 (for d 2 state), respectively. 14 

Figure 1 shows the temperature dependence of the en- 
tropy of the localized state for a series, J p h — 0.01-D ~ 
0.2D with D being the half-bandwidth of the conduc- 
tion band. Other parameters are set as zero in order 
to see an essential character of the model. Then, the 
(16) model Hamiltonian, (14), corresponds to the anisotropic 
Kondo model. For J p h = 0.15D, the entropy of the 



(17) 

U ph EE J ph (I+I- +I-I+), (18) 

H d ee 2(A d -ti + U dc )Il (19) 

The effect of Hu, (6), and Hhyb, (7), will be discussed 
separately below. Here H po t, (16), is a potential scatter- 
ing term. We find that this transformed model Hamilto- 
nian, (14), corresponds to an anisotropic Kondo model 
with a potential scattering and a magnetic (polarized) 
field. Therefore, it is expected that the Kondo effect is 
caused by the pair hopping interaction J p h , if parameters 
are set adequately. 

We study the electronic states of Hamiltonian, (14), on 
the basis of the numerical renormalization group (NRG) 
method of Wilson. 12-14 The model Hamiltonians, (14), 
are transformed into the recursion form as follows: 

Hn+i = A 1//2 i?jv 



+ V, CN+lf N+la fN+la], 



where A is a scale factor (in this letter we set A = 3.0), 
and the coefficients, tjv and ejy, are estimated by the 
tridiagonal procedure. 14, 15 iJjv is given by 

H N = A {N - 1)/2 [U d c + U P h+Ud 

JV-l 

+ T n{fLfn+la + k.C.) 
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Fig. 1. (Color Online) Temperature dependence of the entropy 
for J ph /D = 0.01, 0.1, 0.15, and 0.2, respectively. 



localized state is Si r 



k B lnA at T/D > 10" 2 , and 



(20) g. 



fc s ln2 at 10~ 5 < T/D < 10~ J , and S imp 





at T/D < 10~ 7 . We also find a similar behavior for 
J ph = 0.01D ~ 0.2D. In the region of 5 irap = fc B ln4, 
the localized state is a free state where spin and charge 
degree of freedom are active. On the other hand, in the 
region of Si mp — /csln2, there remains two degrees of 
freedom. In the region of Si mp = 0, the localized state is 
expected to be in the strong coupling singlet state. 

To clarify the nature of the state with S- lmp = fcsln2, 
we estimate the fraction of the occupancy on the lo- 
calized state. Figure 2 shows its temperature depen- 
dence. We use the same parameter set as Fig. 1. Here, 
we focus on the temperature dependence of entropy for 
J P h = 0.15D. At T > 10~ 2 D the fraction of occupancy 
is 0.25 for Q d = 0, 0.5 for Q d = 1, and 0.25 for Q d = 2, 
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Fig. 2. (Color Online) Temperature dependence on the zero-, sin- 
gle, and double occupancy for J p h/D = 0.01, 0.1, 0.15, and 0.2, 
respectively. Q s is the fraction of occupancy. We set the same 
parameters as Fig. 1. 



respectively. Then, we find that the electronic state of 
the localized state is the free state with 4 degrees of free- 
dom. At T < 10~ 3 .D, one finds that the fraction of the 

occupancy of the d° state and the d 2 state are 0.5, while 
the fraction occupancy of the d 1 state is zero. It means 
that the states of d° and d 2 states are degenerate. Thus, 
we conclude that S lmp = fcsln2 stems from degenerate 
charge degrees of freedom (d° and d 2 ), implying that the 
valence skipping state is realized. The same temperature 
dependence is observed also for J p h = 0.01D ~ 0.2D. 

At low temperatures, we find that the energy level 
scheme of the fixed point is exactly the same as that of 
the strong coupling fixed point, although we do not show 
it in this letter. Namely, we find that the ground state is 
the KY singlet state of the charge degree of freedom; the 
charge KY singlet state. 

When Hu and Hhyb are neglected, the model Hamil- 
tonian of eq. (I) or eq.(14) takes the same structure as 
the anisotropic Kondo model with a magnetic field. By 
tuning the chemical potential (/i) and the one-body po- 
tential (Ad) so as to realize T-Ld — 0, the Hamiltonian 
(14) becomes the Kondo model without the magnetic 
field. Then, one finds that the ground state is the charge 
KY singlet by comparing with the above mentioned dis- 
cussions. 

When Tiu and %hyb are considered, the model Hamil- 
tonian of eq. (1) is not the same as the Kondo model. 
However, it is shown that the Kondo effect naturally ap- 
pears by tuning the parameters so as to realize a condi- 
tion Ed« = Edz , where is the energy of the d z state 
(i = 0, 1, 2), estimated by the diagonalization of Hu+Hd 
defined by eq. (6) and eq. (19). E d i (i = 0, 1,2) is esti- 
mated as 

E d o = -A d + fx-U dc , (23) 
E d i = 0, (24) 
E d 2 = U d + A d - fi + U dc . (25) 
Here, we analyze a simple particle-hole symmetric case 



corresponding to — 2Ad — U d . Then, E d <> = E d z is at- 
tained at fx = Udc- We verified that qualitatively the 
same condition is attained, even if we relax the particle- 
hole symmetric condition. 

Figure 3 shows the contour plot of the entropy of the 
localized state in T/D — Vd c /D plane. Here, we choose 
the parameter set as J p h = Udc — 0.1 D, /i = Udc, and 
-2A d = U d = 0.005L>. 




Fig. 3. (Color Online) Contour plot of the entropy of the localized 
state in T/D — Vd c / D plane at J p ^ = Udc = 0.1D, fi = Ud c , and 
-2A = U d = 0.005D. 



The temperature dependence of S lmp for 
Vd c /D < 1CP 2 is almost independent of Vdc/D; 
the electronic state changes from the free state to the 
valence skipping state at T = T* ~ 10~ 2 D, and the 
electronic state changes from the valence skipping state 
to the charge KY singlet state at T = T% ~ 10 _4 Z?, as 
the temperature decreases, where T* is the temperature 
separating the free state and the valence skipping state 
(or the free spin state), and Tr- is the temperature 
separating the valence skipping state (or the free spin 
state) and the charge KY singlet state (the spin KY 
singlet state). For V dc /D > 4 x 1Q _1 , the ground state 
is the spin KY singlet state. As V dc is increased, Tk and 
T* increase. For V dc /D = V^JD ~ 3.659 x 1CT 1 , T K 
drastically decreases, and the entropy is Simp = /csln2 
down to technically zero temperature. From the analysis 
of the energy flow diagram of NRG calculation, we find 
that the electronic state is the coexistence of the charge 
KY singlet and the spin KY singlet. This coexistence is 
possible, because the charge and spin degree of freedoms 
are separated in this model. 

Figure 4 shows the temperature dependence of en- 
tropy, S'imp, and the Sommerfeld coefficient, C- lmp /T, for 
V dc = 0.4D and V dc = 0.3659D. The temperature de- 
pendence of entropy for Vd c — 0.4D and Vd c — 0.36591? 
is consistent with the result in Fig. 3. On the other 
hand, it is found that the Sommerfeld coefficient in- 
creases down to technically zero temperature for the 
seemingly critical value V dc ~ V dc — 0.3659D. It is also 
found that the Sommerfeld coefficient increases in the 
region 10~ 6 D < T < 10" 4 D for V dc = 0.4D, which is 
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Fig. 4. (Color Online) Temperature dependence of entropy, Si m p, 
(black line) and Sommerfeld coefficient Ci mp /T (red line) for 
V dc = OAD and V dc = 0.3659D. 



close to the critical value V5 C . The origin of the enhance- 
ment is suspected to arise from the competition between 
the charge KY singlet and the spin KY singlet. The de- 
tails of the novel fixed point will be discussed elsewhere. 

Finally, we discuss the possibility of the pair hopping 
mechanism in Tl doped PbTe. The density of states 
(DOS) of PbTe is illustrated in Fig. 5(a). 16 PbTe has 
a narrow gap where a direct gap is about 0.17eV, and 
an indirect gap is about 0.025 eV. 17 When a Tl atom is 
doped in PbTe, the DOS of Pbi^Tl^Te is illustrated as 
in Fig. 5(b), where new two peaks appear: one peak near 
the top of the valence band is called a deep state, and the 
other peak much lower than the valence band is called 
a hyper-deep state. 18 The energy level of the deep state 
is the characteristic property in Tl doped PbTe. We find 
that the doping of Tl corresponds to the hole doping in 
the valence band. Indeed, the direct hole doping into the 
valence band has been observed by the photoemission 
spectroscopy. 19 



Valence band 



Conduction band 



l Pbj ,TI,Te 
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IH 


Hole doping 



Hyper- deep d state 
state 



Fig. 5. (Color Online) Schematic picture of the density of states 
in (a) PbTe and (b) Tl doped PbTe. The deep state is located 
at a lower energy level than the top of the valence band. 



It is also known that the hyper-deep state and the 
deep state are almost localized states consisting of the s 



orbital of Tl and p orbitals of Te around the s orbital of 
Tl. 16,20 Thus, we can regard the deep state as the impu- 
rity orbital as in a d or f orbitals, while we neglect the 
hyper-deep state, because the energy level of the hyper- 
deep state is far below from the Fermi level. 

In Tl doped PbTe, the pair hopping interaction is not 
the interaction between two atomic orbitals, but the in- 
teraction between the valence (conduction) band and the 
deep state. The pair hopping interaction can be expressed 
by the linear combination of the atomic- Coulomb inter- 
actions between two p orbitals in Te. It is expected that 
the amplitude of the pair hopping is of the order of the 
parameters as used in Fig. 3. Thus, by tuning the chemi- 
cal potential, it is expected that the charge Kondo effect 
is caused in a realistic parameter set as discussed in Fig. 
3. 

In conclusion, a new microscopic model of the charge 
Kondo effect has been proposed, and its model has been 
studied on the basis of the numerical renormalization 
group method. It has been shown that the charge Kondo 
effect and the valence skipping phenomenon are caused 
by the pair hopping interaction between the conduction 
band and the localized band. It has been also found that 
the charge and spin KY singlet states coexist in this 
effective model for the critical parameter set, and the 
Sommerfeld coefficient is enhanced by the competition 
between these singlet states. 
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